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I Recently, in a series of works by Liu and Li (L&L), they claimed that there exists a timing 

asynchrony of —25.6 ms between the spacecraft attitude and radiometer output timestamps in the 
' original raw WMAP time-ordered data (TOD). L&L reprocessed the WMAP data while the afore- 

mentioned timing asynchrony has been corrected, and they obtained an alternative CMB map in 
which the quadrupole dropped to nearly zero. In the present work, we try to see the implications 
to dark energy cosmology if L&L are right. While L&L claimed that there is a bug in the WMAP 
pipeline which leads to significantly different cosmological parameters, an interesting question nat- 
, urally arises, namely, how robust is the current dark energy cosmology with respect to systematic 

errors and bugs? So, in this work, we adopt the alternative CMB data of L&L as a strawman to 
study the robustness of dark energy predictions. 
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I. INTRODUCTION 

Dark energy cosmology has been one of the most active fields in astronomy and physics, since the 
exciting discovery of current accelerated expansion of our universe [l| . The first evidence came from the 
observation of Type la supernovae (SNIa) in 1998 Five years later, cosmology entered the so-called 
"precision era" in 2003 when the first year Wilkinson Microwave Anisotropy Probe (WMAP) observations 
of the cosmic microwave background (CMB) had been released Up to now, the CMB observation 
is still a very powerful probe for cosmology, and the CMB data from WMAP mission provide the most 
important basis. 

However, in the passed years, some unusual phenomena have been found from the CMB data released 
by WMAP team. Remarkably, among these unusual phenomena, it is claimed that there exists a 
preferred direction in the CMB temperature map (known as the "Axis of Evil" in the literature) [H . In 
fact, there are two different approaches to deal with this problem. The first one is to admit that this 
phenomenon is an observational fact and try to explain it in the cosmological theories (see e.g. 043 
and references therein). The second approach is to consider instead that this phenomenon might be 
an artifact due to the observational systematics. For example, in a series of works by Liu and Li 
(L&L) ^8i-ll0l|. they claimed that there exists a timing asynchrony of —25.6 ms between the spacecraft 
attitude and radiometer output timestamps in the original raw WMAP time-ordered data (TOD). If 
one fixed this problem, the most of CMB quadrupole component disappear. It is worth noting that 
recently their findings has been confirmed by several independent authors (see e.g. [11 13]). In L&L 
reprocessed the WMAP data while the aforementioned timing asynchrony has been corrected, and they 
obtained an alternative CMB map in which the quadrupole dropped to nearly zero. In addition, using the 
alternative CMB data, they constrained the cosmological parameters of the standard flat ACDM model, 
and found that these cosmological parameters have been changed notably. For convenience, we reproduce 
their main results in Table ID It is easy to see that the fractional matter density flmo ~ ^bo + ^co of 
L&L is considerably larger than the one of WMAP [li,[l|. 



Description 


Symbol 


WMAP \U,15} 


L&L 


[8J 


Hubble constant (km/s/Mpc) 


Ho 


71.9^2-6 


71.0 ± 


:2.7 


Baryon density 




0.0441 ± 0.0030 


0.052 ± 0.0030 


Cold dark matter density 




0.214 ±0.027 


0.270 ± 


0.027 


Dark energy density 


Aao 


0.742 ± 0.030 


0.678 ± 


0.030 


Flue. Ampl. at 8/i"^ Mpc 


erg 


0.796 ± 0.036 


0.921 ± 


0.036 


Scalar spectral index 


ris 


q63+° °^* 
u.yoo_Q Qj^g 


0.957 ± 


0.015 


Reionization optical depth 


T 


0.087 ±0.017 


0.109 ± 


0.017 



TABLE I: The cosmological constraints on the standard flat ACDM model, reproduced from Q|. 

To our knowledge, there is still controversy about the findings of L&L in the community. In the present 
work, we try to be neutral as much as possible. Instead of discussing the detailed data process of WMAP 
mission (as in the works by L&L [Si-ilOj or other authors [ll| - [l3| ). in this work we would like to see the 
implications to dark energy cosmology if L&L are right. While L&L claimed that there is a bug in the 
WMAP pipeline which leads to significantly different cosmological parameters, an interesting question 
naturally arises, namely, how robust is the current dark energy cosmology with respect to systematic 
errors and bugs? So, in this work, we adopt the alternative CMB data of L&L as a strawman to study 
the robustness of dark energy predictions. 

As is well known, using the full CMB data to perform a global fitting consumes a large amount of 
computation time and power. As a good alternative, one can instead use the shift parameter R from 
the CMB data, which has been considered extensively in the literature. It is argued that the shift 
parameter R is model- independent, and it contains the main information of the CMB data [1^1 • So, in 
Sec. |TT]we firstly derive the corresponding shift parameter R from the alternative CMB data of L&L Q. 
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In addition, we briefly introduce the other observational data, such as SNIa and the baryon acoustic 
oscillation (BAO), which are also used in the present work. In Sec. IIIIl we perform the data analysis by 
using the cosmological observations introduced in Sec. |TT1 In particular, we discuss the tension between 
CMB and SNIa in Sec. IIII Al we study the age problem in dark energy models in Sec IIIIBI and the 
cosmological constraints on dark energy models are considered in Sec. IIII CI Finally, the conclusion and 
discussions are given in Sec. IIVI 



II. OBSERVATIONAL DATA 



A. Shift parameter R from the alternative CMB data 

As mentioned above, using the shift parameter R from the CMB data is a good approach to fit the 



cosmological models. As is well known, the shift parameter R of the CMB data is defined by [17[ (see 
also ^M^) 



1/2 Z"^* dz 



where flmo is the present fractional energy density of pressureless matter; E = H/Hq in which H = a/a 
is the Hubble parameter; a = (l + z)"-^ is the scale factor (we have set ao = 1; the subscript "0" indicates 
the present value of corresponding quantity; z is the redshift); a dot denotes the derivative with respect 
to cosmic time t; the redshift of recombination is given by [l^ (see also ITsjl) 



= 1048 



1 + 0.00124 (r!fco/i') "'^^^ 1 + 31 (r!„,o/i')^' , (2) 



where f2mo — ^60 + f^cOj ^^nd fihOi ^co are the present fractional energy densities of baryon and cold dark 
matter, respectively; h is the Hubble constant Rq in units of lOOkm/s/Mpc; and 

_ 0.0783 (r>fco/>?)'° ''' _ 0.560 

" 1 + 39.5 (f7,o/i2)" ™' ' " 1 + 21.1 (f7,o/i2)i «i ■ ^ ^ 

The shift parameter R relates the angular diameter distance to the last scattering surface, the comoving 
size of the sound horizon at and the angular scale of the first acoustic peak in CMB power spectrum 
of temperature fluctuations 0, [l3| • 

In principle, one should obtain and R from the full (alternative) CMB data. However, this is a 
hard work consuming a large amount of computation time and power. On the other hand, L&L have 
not published their full alternative CMB data in which the timing asynchrony of —25.6 ms between 
the spacecraft attitude and radiometer output timestamps in the original raw WMAP time-ordered 
data (TOD) has been corrected. Therefore, in this work we use instead the Monte Carlo method (T9l |. 
We choose the standard flat ACDM model to be the fiducial model, since its cosmological parameters 
have been constrained by both WMAP [l3, EH and L&L Q (see Table |l|. We generate the Gaussian 
distributions for rihO; ^co and h from their best- fit parameters and the corresponding Icr uncertainties 
given in Table ID Then, we randomly sample the parameters ilbOj ^^co and h from their corresponding 
Gaussian distribution for iVmc times. For each {r2f,0; ^^cO: we can obtain the corresponding 2* and 
R from Eqs. ^ and ([T]), respectively. Finally, we can determine the means and the corresponding la 
uncertainties for and R from these A'mc samples, respectively. Notice that in this work we have done 
A'mc = 10^ samplings. At first, we check this method by calculating z» and R from the WMAP data, and 
find that our results are well consistent with the ones given by WMAP team [l4|, [Tsl . Then, we turn to 
the alternative CMB data of L&L [8] , and finally find that the corresponding redshift of recombination 
z* is given by 

z* = 1088.610 ±2.373, (4) 

and the shift parameter R is given by 

i?= 1.761 ±0.014. (5) 
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Obviously, the shift parameter R from L&L's alternative CMB data is fairly larger than the one of 
WMAP 111 [H. It is argued that the shift parameter R is model- independent, and it contains the 
main information of the CMB data Therefore, in the literature the shift parameter R has been used 
extensively to constrain cosmological models. For the alternative CMB data of L&L Q , the corresponding 
from the shift parameter R reads 

2 _ f R-1.76l Y 
XcMB- I 0.014 j ■ 



B. SNIa and BAO data 



In addition to the CMB data, we also consider the observations of SNIa and BAO. The data points of 
the latest 557 Union2 SNIa compiled in ^3 are given in terms of the distance modulus Hobsizi). On the 
other hand, the theoretical distance modulus is defined as 

fithizi) = 5log^QDL{zi) + , (7) 

where /xq = 42.38 — 51ogio h and 

in which p denotes the model parameters. The corresponding from the 557 Union2 SNIa reads 

2 / \ [f^obsjzi) - fJ.th(Zi)f , , 

where a is the corresponding la error. The parameter /io is a nuisance parameter but it is independent of 
the data points. One can perform an uniform marginalization over fiQ. However, there is an alternative 
way. Following l 2lL [22l . [28l | , the minimization with respect to /iq can be made by expanding the of 
Eq. © with respect to /xq as 

xl{p)=A-2tioB + filC, (10) 

where 

[fJ'obsiZi) - llth{Zi\pL0 0, p)]^ 



i(p) = E ■ 



^/ X ^ f^obsjzi) - l^th{zi;m = o,p) c ^ 

Eq. (llOp has a minimum for — B /C at 

x^(p) = i(p)~:^. (11) 

Since mm — mm obviously, we can instead minimize which is independent of /io- It is worth 
noting that the corresponding h can be determined hj = B/C for the best-fit parameters. 

Similar to the case of shift parameter i?, one can use the the distance parameter A of the measurement 
of the BAO peak in the distribution of SDSS luminous red galaxies ,23] , which is also mo del- independent 
and contains the main information of the observation of BAO. The distance parameter A is defined by [1^ 



Zb Jo 
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where Zb = 0.35. In [H], the value of A has been determined to be 0.469 (n^/O.gS)-^-^^ ± 0.017. Here 
the scalar spectral index Ug is taken to be 0.957 from L&L (see also Table H]). The corresponding 
xIao = - AobsV/al. ^ _ _ _ ^ ^ 

The best-fit model parameters are determined by minimizing the corresponding ■ As in |25l . l26l |76[ , 
the 68.3% confidence level (C.L.) is determined by Ax^ = — Xmin ^ 2.3 and 3.53 for Up = 1, 
2 and 3, respectively, where Up is the number of free model parameters. Similarly, the 95.4% C.L. is 
determined by Ax^ =x^— Xmin — 4-0, 6.17 and 8.02 for Up — 1, 2 and 3, respectively. 



Observation 






wo 


■Wa 


SNIa 


541.430 


0.420 


-0.863 


-5.490 


SNIa+A 


542.636 


0.278 


-1.007 


-0.105 


SNIa+i? 


542.523 


0.295 


-1.021 


-0.314 


Sm-A+A+R 


542.936 


0.288 


-0.969 


-0.529 



TABLE IL The Xmin and the best-fit values of Q,mO, wo and Wa for the various observations. 




-L75-L5-L25 -1 -0.75-0.5-0.25 _i.4 _i.2 -1 -0.8 -0.6 -0.4 

Wo Wo 



FIG. 1: The 68.3% and 95.4% C.L. contours in the wq — Wa plane and the wq — fimo plane for the observations 
of SNIa alone (blue dashed lines) and SNIa-|-A-|-_R (black solid lines). We also show the best- fit values for the 
observations of SNIa alone (blue box), SNIa-|-A (red star), SNIa-|-_R (green triangle) and SNIa-|-A-|-_R (black point). 

III. DATA ANALYSIS 

In this section, we perform the data analysis by using the cosmological observations mentioned in the 
previous section. In particular, we discuss the tension between CMB and SNIa in Sec. lIII Al we study the 
age problem in dark energy models in Sec IIII Bl and the cosmological constraints on dark energy models 
are considered in Sec. IIII CI 

A. Tension between CMB and SNIa 

There is a long history of the tension between CMB and SNIa. In [27l.l28t , the Gold04 SNIa dataset [29| 
has been shown to be in 2a tension with the SNLS SNIa dataset [30| and the WMAP observations. 
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Although the Gold04 sample was updated to Gold06 [3l| several years later, the tension stih persisted. 



In [32|, the tension and systematics in the Gold06 SNIa dataset has been investigated in great detail. 
Later, the number of SNIa has been significantly increased to 300 ~ 400 when the Union [s^ and 
Constitution [s^l SNIa datasets have been released. Subsequently, in [11] it is found that there is still 
significant tension between the WMAP CMB data and the Union/Constitution SNIa datasets. The Jatest 
SNIa dataset is the Union2 compiled in [2(|, which is an updated version of Union. However, in .3(1 ISTj 
it is found that the tension between the WMAP CMB data and the Union2 SNIa dataset still persisted. 

Obviously, the tension comes from both sides: CMB and SNIa. In (sl. [35l - l37l |. the authors mainly 
concentrated in SNIa datasets, and tried to find the outliers responsible for the tension. However, from 
the works by L&L we become aware of the possible systematics in the CMB data. So, it is natural 

to see whether the alternative CMB data of LfcL jSj can alleviate the tension. 

In this subsection, we follow the discussions in [3a] , and consider the familiar Chevallier-Polarski-Linder 
(CPL) model [38], in which the equation-of-state parameter (EoS) of dark energy is parameterized as 

twde = Wo + 1^0(1 - a) = Wo + iUq— - — , (13) 

1 + z 



where wo and Wa are constants. As is well known, the corresponding E(z) is given by [l9l l39l - l4ll . |71| 

-,1/2 



,0(1 + Zf + (1 - ano) (1 + z)3(l+-o + ^„) exp (^-^) 



(14) 



Here, we consider the shift parameter R from the alternative CMB data of L&L [s*], the latest Union2 
SNIa dataset [23], and the distance parameter A [2^] with from L&L [^. We fit the CPL model to 
the observations of SNIa alone, SNIa+A, SNIa+i?, and SNIa+A+i?, respectively. The best-fit values 
are presented in Table HIl In Fig. [1] we also present the 68.3% and 95.4% C.L. contours in the wq — Wa 
plane and the wq — ^mo plane. From Fig. [U we find that there is still tension (beyond 2a) between the 
alternative CMB data of L&L [8] and the latest Union2 SNIa dataset [13] ■ However, when we compare 
Fig. [Ijand Table HIl of the present work with Figs. 1, 2 and Tables I, II of [11], it is easy to see that the 
tension has been alleviated to some extent in fact. The alleviation mainly comes from two sides: the 
best-fit ri„io for L&L's alternative CMB data is increased, whereas the best-fit J7mO for the Union2 SNIa 
dataset is decreased. It is worth noting that the contribution to the alleviation from Union2 SNIa dataset 
is larger than the one from L&L's alternative CMB data [^. 



B. Age problem 

In history, the age problem played an important role in the cosmology for many times (see e.g. [42| for 
a brief review). The main idea is very simple: the universe cannot be younger than its constituents. For 
example, the matter-dominated Friedmann-Robertson- Walker (FRW) universe can be ruled out because 
its age is smaller than the ages inferred from old globular clusters. The age problem becomes even more 
serious when we consider the age of the universe at high redshift (rather than today, z = 0). So far, 
there are some old high redshift objects (OHROs) been discovered. For instance, the 3.5 Gyr old ga laxy 
LBDS 53W091 at redshift z ^ 1.55 [H, the 4.0 Gyr old galaxy LBDS 53W069 at redshift z = 1.43 [3, 
the 4.0 Gyr old radio galaxy 3C 65 at z = 1.175 [H], and the high redshift quasar B1422+231 at 
z = 3.62 whose best-fit age is 1.5 Gyr with a lower bound of 1.3 Gyr [i^]. In addition, the old quasar 
APM 08279-^5255 at z = 3.91 is also used extensively, whose age is estimated to be 2.0—3.0 Gyr (STj-lisj. 
To assure the robustness of our analysis, we use the most conservative lower age estimate 2.0 Gyr for 
the old quasar APM 08279+5255 at z = 3.91 [13, SI], and the lower age estimate 1.3 Gyr for the high 
redshift quasar B1422+231 at z = 3.62 jig . In the literature, there are many works on the age problem 
in the dark energy models, see e.g. (42l. |49|- 55| and references therein. 

In this subsection, we would like to consider the age problems in the flat ACDM model and the 
holographic dark energy (HDE) model, which have been discussed in the literature by using the earlier 
observational data (see e.g. [42], HO, jEHlll]). Of course, our main goal is to see whether the age problem 
can be alleviated with the alternative CMB data of L&L M . 
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Description 




5(3.91) 


5(1.43) 


5(1.55) 


5(1.175) 


5(3.62) 


best fit 


0.288 


0.799 


1.120 


1.195 


1.310 


1.345 


la lower bound 


0.275 


0.817 


1.144 


1.221 


1.337 


1.376 


la upper bound 


0.302 


0.781 


1.097 


1.170 


1.283 


1.315 


2(7 lower bound 


0.263 


0.836 


1.169 


1.248 


1.366 


1.408 


2a upper bound 


0.316 


0.763 


1.074 


1.145 


1.257 


1.286 



TABLE III: The ratio S{z) = T.^{z)/Tobj aX z = 3.91, 1.43, 1.55, 1.175 and 3.62, for various model parameters 
f^mo (witflin 2a uncertainty) of tlie ACDM model (the corresponding h = 0.696). 



1. Age problem vn the ACDM model 



The age of our universe at redshift z is given by [42, l50|, |5 



It is convenient to introduce the so-called dimensionless age parameter [sO, HJ] 

Uz) ^ Hotiz) ^ r ■ (16) 

At any redshift, the age of our universe should be larger than, at least equal to, the age of the OHRO, 
namely 

Tz{z) > Tobj = Hotobj , or equivalently, S{z) = T^{z)/Tobj > 1 , (17) 

where tobj is the age of the OHRO. It is worth noting that from Eq. p^ . Tz{z) is independent of the 
Hubble constant Hq- On the other hand, from Eqs. ([T7]), Tobj is proportional to the Hubble constant Hq. 
The lower Hq, the smaller Tobj is. 
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1 2 3 4 5 



z 

FIG. 3: The dimensionless age parameter Tz as a function of redshift z for various fimo in tlie ACDM model. 



As is well known, for the flat ACDM model, 

E{z) = Vf^„,o(l + z)3 + (l-r!„,o). (18) 

It is easy to obtain the total — xfi + Xcmb + Xbao ^ function of the single model parameter ilmo 
for the ACDM model, by using the combined observations of the shift parameter R from the alternative 
CMB data of L&L |S] , the latest Union2 SNIa dataset [13] , and the distance parameter A with ris 
from L&L We present the corresponding x^ ^md likelihood C oc in Fig. [2] The best fit has 

Xmin = 545.239, whereas the best-fit parameter is flmo = 0-2881o;q}3 (with la uncertainty) to'o26 (with 
2(7 uncertainty). The corresponding h — 0.696 for the best fit. It is worth noting that although the 
best-fit is smaller than the one from the alternative CMB data of L&L alone [^ (see also Table H]) 
due to the influence of SNIa and BAO, it is still larger than the one from WMAP alone [3) [13 (see also 
Table m and the one from WMAP-|-SNIa-|-BAO [Mlil. 

In Table Hill we present the ratio S{z) = Tz{z)/Tobj at z = 3.91, 1.43, 1.55, 1.175 and 3.62, for various 
model parameters (within 2a uncertainty) of the ACDM model. Obviously, Tz{z) > Tobj holds at 
z = 1.43, 1.55, 1.175 and 3.62, whereas T^(z) < Tobj aX z = 3.91. So, the old quasar APM 08279+5255 
at z = 3.91 cannot be accommodated (beyond 2a). The age problem still exists in the ACDM model, 
even when the alternative CMB data of L&L 8] has been taken into account. 

Let us have a closer observation. As mentioned above, the alternative CMB data of L&L [8] favor a 
larger flmo (this point has also been mentioned by L&L themselves Q). So, it is natural to see how the 
model parameter affects the age of our universe at redshift z. In Fig. [31 we present the dimensionless 
age parameter Tz as a function of redshift z for various ilmo in the ACDM model (notice that from 
Eq. nil), Tz{z) is independent of the Hubble constant Hq). It is easy to see that at any redshift z, the 
larger flmo, the smaller age is. Therefore, the age problem becomes even worse in the ACDM model, 
since the alternative CMB data of L&L [8] favor a larger flmo- 



2. Age problem in the HDE model 

Previously, in [5^, the age problem in the HDE model has been discussed in great detail. Here we 
would like to see whether the age problem can be alleviated with the alternative CMB data of L&L ^ . 
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0.75 
0.7 

0.65 
0.6 

0.55 




0.26 



0.28 



0.3 



0.32 



0.34 



FIG. 4: The 68.3% and 95.4% C.L. contours in the QmO 
parameters are also indicated by a black solid point. 



c parameter space for the HDE model. The best-fit 



The HDE model [56|, |57| was proposed from the holographic principle 58] in the string theory. HDE 
is now an interesting candidate of dark energy, which has been studied extensively in the literature. The 
energy density of HDE reads ^ (see also e.g. (si. [sol, |60| ) 



PA = "ic^m^pL ^ , 

where the numerical constant 3c^ is introduced for convenience, and mp 
Planck mass. In [stI, L has been chosen to be the future event horizon 



dt 

Rh = a J — 



a 



da 



(19) 

BttG)^^/^ is the reduced 
(20) 



From Eqs. 



I ,. „^J1) . and the energy conservation equation + 3Hp\{l + w\) ~ 0, it is easy to find 
that (see e.g~[5l[51 M, Isl) 



dz 



= -(i + z)-ir!A(i-OA) (i + 



(21) 



where f^A is the fractional energy density of HDE. From the Friedmann equation = {p„i + pa) / (S^^ip) 
we have 



1 - ^Aiz) 



1/2 



(22) 



There are 2 independent model parameters, namely ilmo and c. One can obtain fiA{z) by solving the 
differential equation (j2ip with the initial condition ^Aiz = 0) = 1 — f^mo- Substituting i^Aiz) into 
Eq. (22), we can find the corresponding i?(2;) and then the total = + Xcmb +Xbao- By minimizing 
the total x^, we find the best-fit parameters flmo = 0.298 and c — 0.658, while Xmm — 546.628. The 
corresponding h = 0.703 for the best fit. In Fig. |4l we present the corresponding 68.3% and 95.4% C.L. 
contours in the fimo — c parameter space for the HDE model. Comparing with the results of [H^ |63| , it 
is easy to see that flmo becomes larger, whereas c becomes smaller. Note that c < 1 within 2a region, 
the EoS of HDE can be less than -1 il. 
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Description 




5(3.91) 


S'(1.43) 


5(1.55) 


5(1.175) 


5(3.62) 


best fit 


(0.298, 0.658) 


0.771 


1.079 


1.151 


1.263 


1.298 


la left edge 


(0.279, 0.668) 


0.797 


1.112 


1.187 


1.301 


1.341 


la right edge 


(0.319, 0.642) 


0.746 


1.046 


1.116 


1.226 


1.257 


la bottom edge 


(0.303, 0.589) 


0.766 


1.076 


1.147 


1.261 


1.290 


la top edge 


(0.296, 0.734) 


0.772 


1.076 


1.148 


1.259 


1.300 


2a left edge 


(0.267, 0.678) 


0.814 


1.135 


1.211 


1.327 


1.370 


2a right edge 


(0.333, 0.622) 


0.731 


1.027 


1.095 


1.204 


1.231 


2a bottom edge 


(0.308, 0.548) 


0.761 


1.070 


1.141 


1.255 


1.281 


2(7 top edge 


(0.294, 0.789) 


0.774 


1.075 


1.148 


1.257 


1.302 



TABLE IV: The ratio S{z) = T^{z)/Tobj at z ^ 3.91, 1.43, 1.55, 1.175 and 3.62, for various model parameters 
{QmO, c) (within 2a uncertainty) of the HDE model (the corresponding h — 0.703). 




FIG. 5: The dimensionless age parameter Tz as a function of redshift z for various QmO and c in the HDE model. 



In Table ITVl we present the ratio S{z) = Tz{z)/Tobj a.t z — 3.91, 1.43, 1.55, 1.175 and 3.62, for various 
model parameters {flmo, c) (within 2a uncertainty) of the HDE model. Obviously, Tz{z) > Totj holds at 
z = 1.43, 1.55, 1.175 and 3.62, whereas ^.(z) < Tobj at z = 3.91. So, the old quasar APM 08279+5255 
&t z = 3.91 cannot be accommodated (beyond 2a). The age problem still exists in the HDE model, even 
when the alternative CMB data of L&L [8| has been taken into account. Further, comparing with the 
results in [s^], it is easy to find that actually the age problem becomes even worse in the HDE model. 

Again, we would like to see how the model parameters ^Imo and c affect the age of our universe at 
redshift z. In Fig. [SI we present the dimensionless age parameter as a function of redshift z for various 
r^mo and c in the HDE model (notice that from Eq. Tz(z) is independent of the Hubble constant Hq). 
It is easy to see that at any redshift z, for a fixed c , the larger flmo, the smaller age is; on the other hand, 
for a fixed ri„io, the larger c, the smaller age is. However, from Fig. [5j we also find that the influence 
to the age from 51^0 is significantly stronger than the one from c. As mentioned above, the alternative 
CMB data of L&L favor a larger flmo and a smaller c . Therefore, it is not surprising to find that the 
age problem becomes even worse in the HDE model. 
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FIG. 7: The and likehhood £ oc e ''^ as functions of Qrc for the DGP model. 



C. Cosmological constraints on dark energy models 



In this subsection, we consider the cosmological constraints on various dark energy models, by using 
the combined observations of the shift parameter R from the alternative CMB data of L&L [8] , the latest 
Union2 SNIa dataset [13] , and the distance parameter A [l^] with Ug from L&L Q . It is worth noting that 
the cosmological constraints on ACDM model, CPL model and HDE model have already been obtained 
in Sees. IIII Al and IIIIBI In the foUowings, we further consider the cosmological constraints on the other 
four dark energy models. 
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1. XCDM model 



It is well known that in the spatially flat universe which contains pressureless matter and dark energy 
whose EoS is a constant w^, the corresponding E{z) is given by (see e.g. fsol. l60j) 

Eiz) = ^ano(l + z)3 + (1 - 0„,o)(l + z)3(i+-x) . (23) 

By minimizing the corresponding total — xfi + Xcmb + Xbao ' ^^'^ best-fit parameters flmo = 
0.290 and = -1.063, while xLi„ = 543.45. The corresponding h = 0.700 for the best fit. In Fig. H 
we present the corresponding 68.3% and 95.4% C.L. contours in the flmo — parameter space for the 
XCDM model. Comparing with the corresponding results obtained in [Sa, IfiO] , it is easy to find that Vlmo 
becomes fairly larger, whereas w-^ becomes slightly smaller. 



2. DGP model 



One of the simplest modified gravity models is the well-known Dvali-Gabadadzc-Porrati (DGP) 
braneworld model [1^ [6^, which entails altering the Einstein- Hilbert action by a term arising from 
large extra dimensions. For a list of references on the DGP model, see e.g. [63,[6^ and references therein. 
As is well known, for the spatially flat DGP model (here we only consider the self-accelerating branch), 
E{z) is given by ^^^^ 

E{z) = Vf^mo(l + zf + + , (24) 

where fire is a constant. It is easy to see that E{z = 0) = 1 requires 

f^mO = 1 ~ 2-\/ D,rc ■ (25) 

Therefore, the DGP model has only one independent model parameter ri^c- Notice that < Vlj-c < 1/4 is 
required by < Vlmo < 1- It is easy to obtain the total = x)l + Xcmb + Xbao ^ ^ function of the single 
model parameter Vl^c- In Fig. [71 we plot the corresponding and likelihood £ oc The best fit 

has Xmin — 625.2, whereas the best-fit parameter is Vt^c — 0.122lg;Qg5 (with la uncertainty) to^g^g (with 
2(7 uncertainty). The corresponding h = 0.678 for the best fit. From Eq. (1^ r^^o = 0.302 corresponds 
to the best-fit Vt^c- Comparing with the corresponding results obtained in [59l l60j . it is easy to find that 
Q.rc becomes fairly smaller, and hence firno becomes fairly larger. 



3. New agegraphic dark energy model 



The so-called new agegraphic dark energy (NADE) model was proposed in [63, [fij], based on the 
Karolyhazy uncertainty relation which arises from quantum mechanics together with general relativity. 
In fact, the NADE model is an upgraded version of the agegraphic dark energy (ADE) model [68l - [70j . In 
the NADE model, the dark energy density is given by ^6^ ,67] 



Pq 



where n is a constant of order unity; rj is the conformal time which is defined by 

dt f da 
a 



v= - 



(26) 



(27) 



From the Friedmann equation = {pm + Pq) / (Sm^) , the energy conservation equation pm + 'iH pm = 0, 
Eqs. and ((TT)) . it is easy to find that the equation of motion for the fractional energy density of 
NADE, is given by [H [H 



dQ.q 

dz 



-^q (1 - ^q 



3(1 -t-z)- 



(28) 
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FIG. 8: The and likelihood £ oc e"^^/^ as functions of n for the NADE model. 



On the other hand, from the energy conservation equation pq - 
can find that the EoS of NADE is given by (gI, H^I 



3H{pg+pg) = 0, Eqs. (IMl) and (EH), one 



Pq 



3n a 



(29) 



From Eqs. ([^5]). (I?f)) and (P^. we have a very important result, namely, fig = -n^a? /A in the matter- 
dominated epoch (we strongly refer to [1^ for detailed arguments). Thanks to this special analytic 
feature Vlq — r?c? jA — v?{\ + 2:)~^/4 in the matter-dominated epoch, NADE is a single-parameter 
model in practice. If n is given, we can obtain f2g(z) from Eq. (j28|) with the initial condition ^qi^zini) = 
n?{l + Zini)~'^ /A: at any Zini which is deep enough into the matter-dominated epoch (we choose Zini = 2000 
as in [g^), instead of rig(z = 0) = 1 — fimo at z = 0. Then, all other physical quantities, such as 
^miz) = 1 — ^qiz) and Wq{z) in Eq. (|29p . can be obtained correspondingly. So, i^mo = ^miz — 0), 
fiqo = ^q{z = 0) and Wqo = Wq{z — 0) are not independent model parameters. The only free model 
parameter is n in the NADE model. From the Friedmann equation iJ^ = [pm + Pq) / (Sm^), we have 



E{z) = 



Vt,,,a{l + zf 

l-^q{z) 



1/2 



If the single model parameter n is given, we can obtain ^q{z) from Eq. 
^q{zini) = ?^^(l + 2:mi)"^/4 at Zini- Thus, wc get fi„iO = 1 - ^q{z = 0) 



Therefore, we can find the corresponding total x 



(30) 

with the initial condition 
Then, E{z) is on hand. 



XcMB + Xbao- In Fig. [H we plot the total 



X 



^ and likelihood £ oc e as functions of n. The best fit has 



X^ - 



595.001, whereas the best-fit 



parameter is n = 2.693_gQ72 (with la uncertainty) I0143 (with 2(7 uncertainty). The corresponding 
h = 0.681 for the best fit. On the other hand, we find that i7,„o = 0.299, fi^o = 0.701 and w 
correspond to the best-fit n. Comparing with the corresponding results obtained in [5~ 
that n becomes fairly smaller, and correspondingly becomes fairly larger. 



-0.793 



III, we find 



4- Ricci dark energy model 

The so-called Ricci dark energy (RDE) model was proposed in [t^, which can be regarded as a variant of 
the HDE model mentioned above, while its corresponding cut-off L in Eq. ([T^ is chosen to be proportional 



14 



0.36 



0.34 



0.32 



0.3 




0.28 



0.32 



0.34 



0.36 



0.38 



0.4 



FIG. 9: The 68.3% and 95.4% C.L. contours in the QraO 
parameters are also indicated by a black solid point. 



a parameter space for the RDE model. The best-fit 



to the Ricci scalar curvature radius. In [72!|, there is no physical motivation to this proposal for L in fact. 
Later, in [73,] it is found that the Jeans length Rcc which is determined by Rqq = H + 2H^ gives the 
causal connection scale of perturbations in the flat universe. Since the Ricci scalar is also proportional 
to H + 2H^ in the flat universe, the physical motivation for RDE has been found in 'ts'] actually. In the 
RDE model, the corresponding pA is given by (t^] 



PA = 3aml i^H + 2H^ 

where a is a positive constant (when L is chosen to be Rcc in Eq- (|19|), o ne can see that a 
Substituting Eq. (1311) into Friedmann equation, it is easy to find that [s^, [63, [zSl 



(1 + ^) 



4-2/q 



1/2 



(31) 

in fact). 



(32) 



There are two independent model parameters, namely flmo and a . By minimizing the corresponding 
total — xfi + XcMB + XbaO' ^^'^ ^^'^ best-fit parameters Qmo ~ 0.355 and a ~ 0.317, while 
Xmin — 588.303. In Fig. [SI we present the corresponding 68.3% and 95.4% C.L. contours in the f^mo — cx 
parameter space for the RDE model. Comparing with the corresponding results obtained in (sol. l60l|. it 
is easy to find that ftmo becomes fairly larger, whereas a becomes significantly smaller. 



Comparison of dark energy models 



In Table|Vl we summarize all the 7 dark energy models considered in this work. As mentioned above, it 
is common that L&L's alternative CMB data Q favor a larger in all these dark energy models. Here 
we briefly consider the comparison of these models. Following [63|, we adopt three criterions used 
extensively in the literature, i.e., Xmin/^^f ^ Bayesian Information Criterion (BIC) and Akaike Information 
Criterion (AIC). Note that the degree of freedom dof = N — k, whereas N and k are the number of data 



points and the number of free model parameters, respectively. The BIC is defined by 74 1 



BIC 



-2 In A, 



fclnA^ . 



(33) 
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Model 


ACDM 


17" TV T 

XCDM 


CPL 


DGP NADE 


HDE 


RDE 


Best fit 


^mo 0.288 


J2mO = 0.290 


^mo = 0.288 


D,rc = 0.122 n = 2.693 


Qmo = 0.298 


^mo = 0.355 






= -1.063 


wo = -0.969 




c = 0.658 


a = 0.317 








wa = -0.529 










545.239 


543.45 


542.936 


625.2 595.001 


546.628 


588.303 


k 


1 


2 


3 


1 1 


2 


2 




0.977 


0.976 


0.977 


1.120 1.066 


0.981 


1.056 


ABIC 





4.537 


10.349 


79.961 49.762 


7.715 


49.390 


AAIC 





0.211 


1.697 


79.961 49.762 


3.389 


45.064 


Rank 


1 


2 


3 ~ 4 


7 5 ~ 6 


3 ~ 4 


5 ~ 6 



TABLE V: Summarizing all the 7 dark energy models considered in this work. 



where Cmax is the maximum likelihood. In the Gaussian cases, Xmin ~ ^'^^^^max- So, the difference in 
BIG between two models is given by ABIG = Axl,in + AfclniV. The AIG is defined by (TSj 

AlC = -2\nC„,ax + 2k. (34) 

The difference in AIG between two models is given by AAIG = Ax^in + 2Afc. In Table IVl we present 
X^in/dof, ABIG and AAIG for all the 7 models considered in this work. Notice that AGDM has been 
chosen to be the fiducial model when we calculate ABIG and AAIG. From Table fVl it is easy to see that 
the rank of models is coincident in all the 3 criterions (Xmm/'^o/, BIG and AIG). The AGDM model 
is still the best, whereas the DGP model is still the worst. Gomparing with the corresponding results 
obtained in [1^, [lOl, we find that the rank of these models has only slight change. In other words, the 
infiuence to the rank from L&L's alternative GMB data Q is fairly minor. 

IV. CONCLUSION AND DISCUSSIONS 

Recently, in a series of works by L&L [sl-fToj. they claimed that there exists a timing asynchrony of 
—25.6 ms between the spacecraft attitude and radiometer output timestamps in the original raw WMAP 
time-ordered data (TOD). In 0, L&L reprocessed the WMAP data while the aforementioned timing 
asynchrony has been corrected, and they obtained an alternative GMB map in which the quadrupole 
dropped to nearly zero. In the present work, we try to see the implications to dark energy cosmology if 
L&L are right. While L&L claimed that there is a bug in the WMAP pipeline which leads to significantly 
different cosmological parameters, an interesting question naturally arises, namely, how robust is the 
current dark energy cosmology with respect to systematic errors and bugs? So, in this work, we adopt 
the alternative GMB data of L&L as a strawman to study the robustness of dark energy predictions. 

In this work, we found that L&L's alternative GMB data Q favor a larger Q.mO in all the dark energy 
models. As a result, we found that the tension between GMB and SNIa can be alleviated to some extent, 
since SNIa dataset usually favors a large fimo- However, the age problem becomes even worse in the 
dark energy models, since a larger fimo usually leads to a smaller age of our universe at any redshift z. 
On the other hand, we found that L&L's alternative GMB data % do not significantly change the rank 
of dark energy models from the perspective of model comparison. Of course, it is a big advantage that 
the quadrupole dropped to nearly zero in the alternative GMB map of L&L (see [8l-[To|). Altogether, we 
consider that it is better to keep neutral to L&L's findings so far. While WMAP has been critical in 
establishing a highly successful cosmology, it is most important that these results are verified by other 
independent experiments, such as Planck. 

Finally, we stress that this work is a hypothetical study in fact, since the works of L&L are 
highly controversial in the community. One should be aware of the risk that the present work is likely 
to be completely irrelevant if it turns out that the works of L&L are indeed flawed. Here are some 
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remarks to be taken serious (we thank the anonymous referee for pointing out these issues). Firstly, in 
fact there are many discussions on L&L's claims in the community (see e.g. [77l - l79j ). and many people 
(including both WMAP and Planck experts) believe that L&L are mistaken in their claims. Secondly, if 
the high-£ spectrum really is as discrepant as claimed by L&L, one could expect that many ground-based 
experiments would also have seen this effect. Thirdly, it is worth asking whether the timing offset was 
applied correctly in the actual map making but not in the calibration step. Up to now, to our knowledge, 
the discussions on L&L's claims in the community are still unsettled, whereas L&L are still defending 
their claims persistently in more papers (see [80l | for instance) and in many conferential talks around the 
world. We hope that this ongoing controversy could be finished with a firm and reliable conclusion in 
the near future. 
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